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Abstract 

We construct a frame of complex Gaussians for the space of L 2 (R n ) 
functions. When propagated along bicharacteristics for the wave equa- 
tion, the frame can be used to build a parametrix with suitable error 
terms. When the coefficients of the wave equation have more regularity, 
propagated frame functions become Gaussian beams. 



Introduction 

In A Parametrix Construction for the Wave Equation with C ' Coefficients 
Hart Smith constructed a parametrix solution for the wave equation using a 
frame that is now called curvclcts. In this paper we construct a new frame out 
of Gaussian functions. When a Gaussian function is propagated along the ray, 
it becomes a Gaussian beam which looks like a Gaussian distribution on planes 
perpendicular to a ray in space-time. While the existence of such solutions 
has been known to the pure mathematics community since the 1960's, there 
has been a recent revival of interest given their robustness in approximating 
solutions to PDE's. 

Recently Nicolay Tanushev numerically simulated mountain waves with a 
high degree of accuracy using superpositions of high frequency Gaussian beams. 
Because Gaussian beams are concentrated along a single ray, it is desirable to use 
many of them to represent a solution because a global solution is rarely concen- 
trated along a single curve. Tanushev's thesis showed that Gaussian beams have 
several major advantages over other techniques used to numerically approximate 
the solution to a mountain wave. Motivated by these numerical calculations we 
will show that a frame consisting entirely of complex Gaussians can be used to 
build an accurate parametrix to the wave equation. The idea of using complex 
Gaussians to build a parametrix is not new. In his paper Strichartz Estimates 
for Operators with Nonsmooth Coefficients and the Nonlinear Wave Equation, 
Daniel Tataru constructed a parametrix to the wave equation with low regular- 
ity coefficients using a modified FBI transform. While the solution in his paper 
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is elegant, numerical calculations with such a construction would be difficult if 
not impossible. Representing initial data in terms of a frame of Gaussians may 
lead to more viable and accurate numerical solutions. 

For the rest of this paper we will consider the wave equation 

dfu(t,x) — A(t,x,d x )u(t,x) — d^u(t,x) — aij(t, x)d Xi d Xj u(t 1 x) = 

l<z,_7<n 

and we let 

A(x,t) = {aij(x,t)}i< it j< n . 

We assume that the matrix A is uniformly positive definite and bounded, that 
is there exists a constant C > with 

^< E MM)£^<q£l 2 

l<i,j<n 

for all (t,x,£) in [— T, T] x E" x K™. Here T is fixed and finite. Furthermore 
we assume the entries of the matrix, dij(x, t) with (t, x) G [— T, T] x M n , are in 
C ' . C ' coefficients are of interest because they are minimally regular. They 
satisfy a Lipschitz condition in x and t 

\aij(t,x) - Oij(t',x')\ < C(\t-t'\ + \x - x'\) 

and their first derivatives in x satisfy a Lipschitz condition 

\V x a tj (t,x) - V x a,ij(t,x')\ < C\x - x'\ 

This paper is divided into three major parts. The focus of section 1 is the 
introduction of a frame of Gaussian functions, which will represent elements of 
the Hilbcrt space L 2 (M"). Section 2 introduces Gaussian beams and contains 
the necessary estimates for the construction of a parametrix for the wave equa- 
tion with C ' coefficients. Finally section 3 contains the actual parametrix 
construction. 



1 Construction of the Frame 

Let the set of functions {^ 7 (x)} 7e r be defined as follows: 

7 (x) = ( ^T^ 2 ) ex P («6y ' ( x ~ x i) ~ \^i\\ x ~ x i?) 

where 7 is the index 7 = (i, k, a), i G Ik, where Ik is a finite subset of integers 
that depends on k, k G N, and a G Z™. In the first two lemmas we will pick 
£ 7 = 2 h u>i.k a vector in M. n with 4 < \u)i t k\ < 1 and x 7 = Ax 7 a another vector 
in R™ with Ax 7 a scale factor depending on k. We will show that these vectors 
can be chosen so that the set of functions {</> 7 (a;)} 76 r form a frame for L 2 (R n ). 
Not only will our chosen set of {</> 7 (a;)} 7e r form a frame in L 2 (R n ), but we 
will show that weighted sequences of frame functions are comparable to the m th 
Sobolev norm (provided it exists) of any f{x). In particular: 
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Theorem 1. For any finite m > and f(x) £ H m (M. n ) there exist constants 
C\ and C2 independent of 7 with < C\ < C2 such that the following holds: 

< C X \\f{x)\\% m{VLn) < E I2 fcm c(7)| 2 < C 2 \\f(x)\\% HRn) (1.1) 

7 



c il) = I (t> 7 (x)f(x)dx. 

For this paper we will use the convention that the Fourier Transform for a 
function h(u) £ L 2 (M. n ) is defined as 



Hv) = j e- l7ru h(u) du. 

R" 

We will also need to introduce the following functions 

^l( w )=[^-J ex P (^7 - W ~ l£yllH ) ■ 

The only difference between ■0 7 (x — it) and </> 7 (x) is that the discrete variable 
x 7 is now a continuous one, u, and there is no factor of (Ax 7 )t. Here we note 
that 

|S(f)P = 2-»e x p(-^£). 

In Lemma [T] we construct an approximate partition of unity from the sum of 
the squares of the Fourier transforms of the jjj^(w). 

Lemma 1. One can chose ui^k, i £ Ik, k £ N with | < |wi,fc| < 1 so that the 
inequalities 



< c[\e m < E 22femcx p (- ^ 2 | 2 t^f j - c ^ 2m ( L2 ) 



hold for all £ G M n , m > 0, finite. Here C[ and C ! 2 are constants independent 
oft- 

For clarity, we will save the proof of Lemma [1] and Lemma [2] below until the 
end of the proof of Theorem [TJ For Lemma [2j we will pick Ax 7 so that we can 
approximate the center term in the inequality ([l.ljl by an expression which no 
longer involves a, effectively turning the summation over a into an integral. 

Lemma 2. For fixed k £ N, let Ax 7 be equal C e 2~^~ ek with e > and C e a 
small constant independent of k and dependent on e. Then for every e > there 
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exists a choice of C e such that the following holds 



2 2km f{x)(l) J {x-aAx J )f(x')(j) 1 (x' -aAx 7 )dxdx' 



J2 / 2 2km ^{x-u)f(x)^(x' -u)f{x')dudxdx' 



Proof of Theorem [3 If wc let 



< 



2 2fe "V 7 (x - u)V> 7 (x' - u)f(x)f(x') du dx dx' 

then the kernel of this expression can be rewritten as 

2 2km ^(x - u)^{x' -u)du = (2tt)™ / e l S-( x - x 'h 2km \i^£j\ 2 d£ 



since the Fourier Transform is an isometry on L 2 (M. n ). As remarked earlier, 



|^(0| 2 = 2-"cxp(^- 
so that by Lemma [1] and Fubini's Theorem 

2 



7T 



i^r/(o 



j J 2 2km i>ii% - u)<ip 7 {x' -u)f{x)f(x')dudxdx' 



ier/(e) 

which is equivalent to 

K n C[\\f(x)\\% m{Mn) <J2j f J 2 2km ^(x - u)^(x> - u)f(x)f(x') du dx dx' 

(1.3) 

<^\\f(x)\\% m(Rny 
From Lemma [2] we also have 



53 / / 2 2km f{x)(j) 1 {x-aAx 1 )f{x l )(j) 1 {x l -aAx^dxdx' 



(1.4) 



J2 J I I 2: 2krn >ip 1 (x-u)f{x)il; 1 {x' -u)f(x')dudxdx' 



< 



^"e" 1 „ „ l|2 



L 2 (R" 



but since ^5- < C[ and C 2 > C[ we can combine inequalities (|1.3I) and (|1.4[) to 
conclude 



Ci\\f\\\ mm < 



<C 2 ||/|| 2 ^ 
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2 2km f(x)cj) 1 {x - a; 7 )/(a;')</> 7 (a; , - x 7 ) di' 



which is the result (jl.ip . □ 

Proof of Lemma\7^ Since £ € R™/{0}, we begin by considering R" as an infinite 
union of dyadic annuli, each of which we will cover with real Gaussians which 
are centered at our choice of 2 k uj it k- In every annulus 2 k ~ 1 < |£| < 2 fc , for all 
k S N, we choose the vectors 2 k uji^ such that for all i ^ j \2 k u)i^~ 2 k Uj ! k\ > 2 2 . 
but that the number of 2 k Ui.k in each annulus is as large as possible. The index 
set Ik is finite as the volume of every annulus is finite. 

Fixing £ for the rest of this proof, £ must lie in an annulus 2 fe_1 < |£| < 2 k 
for some fixed k in N. As a result of our choice of vectors, for all £ € R ra there 
exists at least one point 2 k 0Ji } k for which the inequality |£ — 2 k Ui^\ < 2? holds. 
This condition gives the existence of a lower bound: 



I* 



>i,k\ 



iere-<^2-ex P 

(»,fc) 

To show the sum is bounded above, we will consider sets of indices A, B, C, D, 
and E whose union contains all the indices (i, k) in 7 and show that the contri- 
bution to the sum from each of these sets is bounded by a constant multiple of 
|£| 2m . The cases k = 0, 1 are easy, so we consider k > 2 

First let A consist of those indices (i, q) for which |£ — 2 q uji, q \ < 2?. Clearly 
q can only be equal to k — 1, fc, or k + 1. Fixing q for the moment and setting 
r = 2^, if we consider a ball B of radius | centered at each 2 9 Wi jg then for all 
pairs (i,q), (j,q) G A with i ^ j, B(2 q U) i>q , §) n B(2 q uj jtq , §) = 0. But by the 
triangle inequality, all balls of radius | with centers that have indices in A arc 
contained in a ball of radius 4^ around £. Therefore the total number of balls 
N is bounded as 



3r 

2 J ) -' •-" V" V" 2 



Vol ( B ( f , — ) ) > iWoZ ( B ( 2 9 Wi )9 , - 



which implies TV < 3™. Since there are only three possible values q can take, 
the total contribution for the set A to the sum is bounded by 3 ra+1 . 

For the second set let B consist of those indices (i, q) for which the inequality 
2^ < |£ — 2 q u>i_ q j < 2 k holds and also \k — q\ < 1. We can write B as a collection 



of subsets Bj such that 



2~ 

B=\jB j 

J=2 
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where Bj denotes the set of indices for which (j — 1)2 s < |£ — 2 q uj i-q \ < j2 2 . As 
before, we consider balls of radius § = 2s -1 centered at each 2 q uj ijq such that 
for all pairs (i, q), (j, q) e Bj with i % j, B{2 q u] uq , §) n B(2 q u) jtq , §) = 0. By the 
triangle inequality, all balls with centers that have indices in Bj are contained 
in an annulus centered about £ with inner radius (j — l)r — ^ and outer radius 
jr + § . The total number of indices for fixed q in each set Bj is bounded as 

Vol (B (e,jr+|)) - Vol (B - l)r- 0) > NVol (b 0) 
which implies 

«*(H)"-H)> 

Since g can take only three possible values, multiplying this last bound by 3 
gives a bound on the total number of indices in each set Bj. Because of the 
restriction on the size of |£ — 2 q uji^ q \ and the fact that \q — k\ < 1, the inequality 

2 - 2\2 q uj,, q \ ~ 
holds for each tuple in Bj. Summing over all of the sets Bj 

3 (i, 9 )eBj v 1 l - ql / 

<gw)((i4)"-0-i)X- fl ^) 

00 

< ^2 2fem 3(2")(j + l)"exp 

The sum 

00 

^3(2")( J + irexp 

is finite; furthermore it is uniformly bounded regardless of the choice of k and 
hence of £. Therefore since 2 fe_1 < |£| < 2 k , the total contribution from the set 
B is bounded by a constant times |£| 2m . 

Next, let C be the set of indices (i,q) for which |£ — 2^^^ | > 2 k holds 
and also \k — q\ < 1. As before for each fixed q, we take balls of radius | = 
2^ _1 centered at each 2 q iOi^ q so that for all pairs (i,q),(j,q) € C with i 7^ j, 
B(2 q u>i, q , |) C\B(2 q ujj tq: |) = 0. All balls with centers that have indices in C are 
contained in an annulus centered about the origin with inner radius r 2 — ^ and 
outer radius r 2 + Since we have removed a number of the vectors because 
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their indices are in B, the total number of indices, N for fixed q is over-estimated 
as follows 

Vol (B (o,r 2 + 0) - Vol (B (0,r 2 - 0) > WoJ (i? (2V, g , 0) 
which implies 

«-(H)*-H)> 

Since |£ — 2 9 u; l:9 | > 2 fe for all e C, the inequality 

29+ 1 - 2|2« Wl , 9 | 

holds for each point 2 9 w; i? with indices in C. Then the contribution from the 
set C is bounded in terms of a sum over k as 



But since 




and 2 fe ™exp (— 2 fc ~ 2 ) — > for all finite n G R as fc — > oo, (| 1 . 5|) is bounded is 
independent of £. So we can conclude that since 2 k ~ l < |£| < 2 k , the total 
contribution from the set C is bounded by a constant times |£| 2m as well. 

Now let V be the set of indices (i, q) for which q < k — 1. To find the number 
of vectors in V for fixed g we again take balls of radius £ = 2^ _1 centered 
at each 2 q Ui, q so that for all pairs (i,q), {j,q) € V with i ^ j, B(2 q uJi tq , |) (~1 
B(2 q u>j^ q , |) = 0. By the triangle inequality, all balls with centers that have 
indices in T> are contained in an annulus centered about the origin with inner 
radius r 2 — § and outer radius r 2 + § . The total number of indices N is bounded 
as 

Vol (B (V 2 + £)) -Vol (B (0,r 2 - r ~y)>NVol (b (2 9 ^, 9 ,0) 
which gives 

w ^(H)"-H)> 
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We can conclude there are at most 2™ ^ (2 2 + |) " — (2 2 — ^) vectors for fixed 
q. Since for these indices q < k — 1, the inequality 



2 ,-i< (2 fc - 1 -2^) 2 < |e-2^ g | 2 



29+ 1 " 2|29 Wi , 9 | 



holds for each (i, q) in T>. The total contribution from the set V is also bounded 
by a constant times |£| 2m 

£ 2 2fem cxp^ He ~ 2 ^ 



(i, 9 )ex> ^ 1 hql 

< ]rV™2" (( 2 * + ( 2 * - 0") exp(-2^ 1 ) 

CXD 

< ^2 2fcm 2"(2* + l)"exp (-2 9 - 1 ) 
since the sum 

00 

^2"(2* + l)"exp(-2«- 1 ) 

9=1 

is uniformly bounded with respect to k . 

The final set 8 contributing to the sum consists of the indices (i, q) for which 
q > k + 1. Again, as above the total number of vectors N for fixed q is at most 

2" [ [ 2* + ^2* 



2/ V 2 



To find the exponential contribution for each q > k + 1 note that for each 

(»,?)ef 

2? _ 5 < (2^ 1 -2^- 2 ) 2 < |e-2y 9 | 2 _ 



29+ 1 " 2|29w ii9 | 



Therefore 



£ 22fcWex p( 'LCT 12 ) < £ 2 2fe ™2"( 2 f + l) n exp(-2»- 



The sum 



00 

^ 2" + l) exp(-2«- 



g=fc+2 



is convergent and bounded independent of k and q. Therefore the total con- 
tribution from the set 8 is bounded by a constant times |£| 2m as well. This 
completes the proof of the Lemma. The construction of the approximate parti- 
tion of unity is similar in idea to the construction of almost orthogonal frames 
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in Meyer's book Wavelets. The £ 7 which are further away from the variable £ 
contribute less to the the partition than those which are close. 

□ 

Proof of Lemma\M For convenience we let: 

5 7 (w, x, x')(A y x) n = 2 2fc "V 7 (x - u)(j> 7 (x' - u) 
which implies that the operator 

^ ^2 2km cl3-y(x - aAx 1 )4> 1 (x' - aAi 7 ) 

7 

is equal to 

We will rewrite the sum over a above using the Poisson summation formula. 
Recall 

Theorem 2. Poisson Summation Formula: Let a be constant, h(u) G <S(M"), 
and a, /3 € Z" Then the following holds 



a£Z" /3eZ" 

Since by definition 



g^(j),x,x')= / e lr) ' w g 1 (u, x, x') du, 



so we start by computing 
g 7 (u,x,x') 



22km j Ivy 



" ( 2*) ° XP _ ^ ' ^ GXP ^ 2 " 2 + 2U(:E + ^ ~ ^ ~ ^ 
= 2 2 ^ (M)" ex P (i(x-^) .e 7 )exp (-2|£ 7 | (« - (^)) 2 ) CX P (l£rl " ^ + 

= 2^ (M)"exp( l( ,-,0-C 7 )exp(-2|C 7 | ( u - (^)) 2 ) cxp (-M^l! 

which by a standard result on the Fourier transform of a Gaussian (sec Appendix 
1) gives 

g^(r),x,x') (1.6) 



2|£ 7 i; v 2tt y ' v v ' 7 ' V 2 j j * v sie 
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Now we notice that 



g 1 (0,x,x') = / g y (u,x,x')du 



so in applying Poisson summation formula we obtain 

E (g-y(aAx 7 ,x,x'))(Ax 7 ) n = / g 1 {u,x,x')du + ^ h 

aSZ™ J n /3eZ™/3#0 



27T/3 



where g 7 (rj, x, x') is given explicitly by (|1.6[) . From this we can conclude the left 
hand side of (11.41) is 



E E 9j (^-,x,x' ^ f(x)f(x')dxdx' 



By symmetry of the integrands in x and x', if we use Schur's lemma then the 
inequality in (|1.4[) follows from the estimate: 



sup 



E, 27T/3 



(1.7) 



(i,k) P<EZ n ,/3=£0 

If we examine the argument of the integral in (|1 .7[) we find from equality (|1.6j) 



< E 

/3eZ™,/3#0 



•9-; 



27T/3 



(27T/3) 2 I^Kx-x') 1 



8|e 7 |(Ax 7 )2 



Integrating both sides of the above inequality with respect to x' gives 



sup 



/ j Q~t \ a ,x,x 



/3ez™,/3#o 



\ Ax 



exp 



/3(EZ",/3#0 



(27T/3) 2 



8|C 7 l(Aa; 7 ) : 



Let /3 = /32, Pn) then since with this notation /3j € Z is indexed indepen- 
dent of Pj G Z for all i ^ j we have 



e n 



exp 



(2ttA) 2 



8|^l(Ax 7 ) 2 



n e 



exp 



(2ttA) 5 



^ 7 |(Ax 7 )^ 



As /? ^ 0, at least one of the fi^s must also be nonzero. Without loss of generality 
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take (3 n ^ then 

"8|{ 7 |(AxJ 



8|^|(Ax 7 )V;; ^A.*. V 8|{ 7 |(Ax,)^ 
To put a bound on this last expression, we now need to pick Ax 7 . Letting 

m ( 27r ) 2 
aW ^8\^\(A Xi y 

since Aa; 7 is of the form C e 2~^ ek we have that 

7T 2 2 £fe 

»<" - W 

Now if we pick C e < 4 then for all k, a(k) > 1 always. Therefore for any such 
choice of C c we have 

E exp (-iifcO < - t^w < 2 



8|£ 7 |(Az 7 )s 

which ensures the first product in (|1.8p is uniformly bounded independent of /e: 

(27r/?j) z \ \ \ / 2 \ / 2e 



n e«p 

t=l \/3iGZ 



< T TTT < 



|^ 7 |(Aa; 7 )V / / \l-e-°( fc V V 6 " 1 



From Lemma 1, the number of the 2 fc w.; ; fc can be over-estimated by by 2 n (2 2 + 
1)" per fixed k, so combining estimates 



2e x "" ] 



(i,fc)/3eZ™,/3/0 
oo 

<^2 2fcm 7i(2* +1)" ( -^-) cxp(-a(fc)) 
fc=i ^ e ' 

Since e > 0, exp(— a(k)) dominates any power of 2 k . Thus as long as C e is chosen 
sufficiently small we can make this sum less than y ^"^ 1 which concludes the 
proof of inequality (|1.7|) . 

□ 
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2 Operator Norm Estimates 



From Theorem 1 in Section 1, the operator P™(f(y)) = {c(7)} 7 er where 



c( 7 )= / 2* m ^(x)f(x)dx 

is a one-to-one bounded mapping of H m (W 1 ) into the space of sequences weighted 
with 2 km which are convergent in / 2 (r). Let P™ = P™* be defined as follows 

P 2 m : ; 2 (r) L 2 (R n ), P 2 ({c( 7 )}) - ^2 fcm c(7)^(y) 

7 

Now recall that T is an operator of order m if T maps H r (R n ) -> iP" _m (R"). 
In section 1, we showed that 11™ = P™ o P™ is an operator of order 2m. Let 
/ denote the identity operator, as there exist constants C[ and C 2 such that 
C[I < n° < C 2 I, in P 2 (R n ) norm sense, P® is bounded and invertiblc on its 
range. The construction of Pf and P 2 allows us to translate the characterization 
of functions and operators in H m (M. n ) to the framework of weighted sequences 
in l 2 (T). Armed with the frame operators, we will show that when the frame 
functions are propagated along bicharacteristics for the wave equation, their 
Sobolev norm is preserved. This will help us also show that the the action of 
the operator T(x, t, d x , d t ) = df — A(x, t, d x ) on the parametrix is order 1. The 
estimates established in this section will ultimately be useful in building the 
parametrix in section 3. 

First we recall that T(x,t,d x ,dt) = df — A(x,t,d x ) has principal symbol 
p(x, t, £, r) = t 2 — ^2 a,ij(x, The bicharacteristics associated to p are 

dt dx-j d^j dr 



Setting 9 = I X] a i.jiiij i we find that p = (T—q){T+q). There are two choices 
V'i . / 

for null bicharacteristics. Here we assume that r = q 1 so the bicharacteristic 
equations (|2.ip become 

dt dx p< ,„ „, 

S=^ = 2r = 2(? < w = ! = -*' (2 - 2) 

d£ -p x dr 

1x, -tt = 1 



dt 2q dt 



Define 



as the solution to the system (|2.2[) at time t with initial conditions 

(X~f{t ,t , X-y, £7)) &y(t )t 1^7:^7)) = (2-71^7) 
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where (x 7 , £ 7 ) are given in Lemmas Q] and [3] of Section 1. We let U(t, t') denote 
the the evolution operator associated to this transformation. Often we will 
abbreviate 

W(i,0)(a; 7 ,£y) = (x 7 (t,0,£ 7 ,£ 7 ),£ 7 (t,0,x 7 ,£ 7 )) 

as 

(x 7 (t),^(t)) 

and 

li(0,t)(Xy,£y) = {Xy(0,t,Xy,^y),^y(0,t,Xy,^y)) 

as 
Let 

^(t, x) = ( M^J. ^ 2 ex p (if(f) • (x - x 7 (t)) - |f 7 (t)||x - x 7 (t)| 2 ) . 

Then define E(t) to be the propagation operator acting on f(x) <G L 2 (M. n ) as 
follows: 

n°£;(t)n / = p§B E (t)P?f 

= X>E(7,7' ) f)c(7>r0'0 

7,7' 

where 



6^(7)7')*)= / (j) 1 {x)cj)^{t,x)dx 

denotes the entries of the matrix Bsit)- Then note that TE(t) is defined by 
the following equation 

n°T£;(i)n / = p$B T (t)p?f 

= J2b T h,l',t)c( 1 ')^(x) 



where 



br(l>7',t)= <t> 1 (x)T(f)y{t,x)dx 

denotes the entries of the matrix -Bt(^)- The central theorem of this section is 
then: 

Theorem 3. E(t) is a bounded operator of order ; and TE(t) is a bounded 
operator of order 1 . 
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From Section 1, II is bounded and invertible and also by Theorem [T] we 
know the relationship of the frame to the Sobolev norm of f(x). Therefore to 
prove Theorem 3 by Schur's lemma it suffices to show: 

X)iM7,v,t)i<tf, 5> B (7,y,*)i<c (2.3) 

7 7' 

and 

$>t(7,7',0I<^\ £M7,7',*)I <C2 k . (2.4) 

7 7' 

where C denotes a constant independent of 7 and 7'. We will also show this 
constant is uniform for all £ € [— T, T]. 
We start by examining T<pj(t,x): 

Lemma 3. 



V(t, x 7 (£), £ 7 (£)) = £(£) • (a; - z 7 (£)) + i|£ 7 (t)||a: - x 7 (£)| 2 



and 

p(t,a; ) Vx,V't) = 0(|^(*)| 2 | a; - a; 7W| 2 )' 

Proof. As p(t,x,ip x ,ipt) is positive homogeneous of degree two in |£ 7 (£)|, the 
desired conclusion will follow if on nul-bicharacteristics (£, x 7 (£), £ 7 (£)) we can 
show that 

V^(£,x,^(£,x 7 (£),e 7 (£)),^(£,x 7 (£),C 7 (£))) = 0. 
Computing V x p(t,x,i/} x ,ip t ), 
d 

Q^-P(t, x , tpx,1pt) = Px 3 + P&^XtXj + Vri>T X] ■ (2.5) 

Dividing (|2.5|) by 2q and substituting the equations in (|2.2j) into the right hand 
side of (|2.5[) we obtain 

-1 + ^,^.+^ (2.6) 

As t/;^ ■(£, Xj(t), £ 7 (£)) = £j(£), differentiating £j(£) with respect to £ we have 

d£w da;/ , .„ „. 

1= < 2 - 7 > 



Substituting (|2.7[) into (|2.6[) implies (|2.6p is 0, which happens if and only if (|2 .5[) 
vanishes on nul-bicharacteristics. □ 
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With Lemma [3] in mind, we consider the entries of the matrices -Be(£) and 
Brit). First we set 



0° , 



\^\\^{t)\Ax 7 Ax Y 



then 

&b(7:7'j*) 

= /3° j7 , / exp (i(x — x~f'(t)) ■ £y(i) — i(x — x 7 ) ■ £ 7 — |£y(i)||:r — Xj>(t)\ 2 — |£ 7 |||o; — £ 7 | 2 ) dx 



and to leading order 



(2.8) 



/3° 7 , / exp (i(x - Xy>(i)) ■ £ 7 >(t) - i(x - x 7 ) ■ £ 7 - |£ 7 '(f)||x - xy (i)| 2 - |£ 7 ||x ■ 



x 7 | 2 ) 



x |x-x r (i)| 2 |£ r (t)| 2 dx. 

The first inner product, 6^(7, 7', t) is evaluated 
M7>7',*) 

= /3° >7 , / exp (i(x - xy(t)) • ~ ^ ~ ^7) ' £7 ~ 1 1 21 ~ x y( t )\ 2 - \tj\\\ x ■ 



xJ 2 ) dx 



0°,-y> eX P ( ix l ■ £7 ~ ix l' (*) ' £7' (*)) 

exp (ias • - f 7 )) exp f - (|£ 7 | + |f 7 /(t)|) 

lowlier 



|£ 7 |x 7 + |f 7 /(i)|x 7 /(i) 



K7I + M 



x exp 



iSrl + IM*)! 
Making the change of variable 



\xj'(t) — Xj\ I dx 



y = x 



l£ 7 |z 7 + |g 7 '(t)|s 7 /(t) 
ISrl + IM*)! ' 



(2.9) 



6e(7, 7', i) takes the form of the Fourier transform of a Gaussian integral which 
one can evaluate (see Appendix 1) to obtain: 



/3 7 , 7 ' exp ( i ( x 7 ■ £ 7 — x 7 < (t) ■ £ 7 < (t) + 
2 



|£ 7 |x 7 + |f 7 /(i)|x 7 >(i) 



x exp 



4(ie 7 i + ier'(*)i) 



exp 



I£ri + M*)l 



(er'(*)-£r) 



15 



where now 



|g 7 ||g r ft)|Ax 7 Ax r 
4 7 r(|£ r ,(i)| + |£ r |) 



(2.10) 



so that 

\bE{j,l',t)\ 

< /3 7 . 7 ' x exp 



(2.11) 



4(|£ r | + |£r'(*)l) 



exp 



|£ 7 I + IM*)I 



For the integral &t(7, 7',i) we make the same substitution (|2.9j) into f|2 . 8[) . 
Then we set 

»? = fy(*)-^7. c=|6rl + M*)l 

and 

|g 7 |(x 7 - ay ft)) 
l£rl + l£r'(*)l 

so that we can apply the estimates in Appendix 1 which gives [brij,^ ,t)\ is 
equal 



/^'M'exp ( ~ 



rj 2 ibr\ 
Ac 2 c 



— + b 2 + — 



1 

2c 



Applying Cauchy- Schwartz and back substituting values for 77, c, and 6 so that 
for C a constant independent of 7,7' 



M7,7',*)l <<?/3 7 ,y 



IM*)I%-M*)I 2 1 Mt)| 2 |$ 7 | 2 |x 7 -Mt)| 2 
(IMOI + K7I) 2 (N + MOI) 2 



x exp 



IMt)-g 7 | 2 

4(U + |^(*)|) 



exp 



lerl + IM*)l 



< 



C/3 7 , 7 ' + IM*) 1 1^7 1 1*7 - ^7'WI 2 ) 



x exp 



4(|£rl + |M*)l) 



exp 



|£ 7 I + M*)I 



jiij'-y' X-y I 



The next two Lemmas characterize properties of the evolution operator 
U(t,t') acting on the lattice which will assist us in obtaining the bounds 
and 



Lemma 4. As in the introduction we let 

A(x,t) = {aij(x,t)}i<i t j<„ 

be a real symmetric n x n matrix with entries a,ij(x,t) in C ' . For the rest 
of this lemma C > denotes a constant which is independent of the essential 
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variables. Furthermore, as before A(t, x) is bounded and positive definite and 
we let 

If we consider the system 



• ■j 



dx 



(2.12) 



with initial conditions \x(0)\ < R and - < |£(0)| < a for some finite a, R > 
then solutions to the system h2.12\) satisfy the following two conditions. 

1. |a?(i)-a?(0)| < Cy/n~\T\ 

2. For all finite T > there exists a constant C(T,a) such that 

< \u(t)\<C(T,a) 



C(T,a) 
whenever \t\ < T 
Proof. We prove condition (1) first and then condition (2). 
1. Computing q^ 



d_ 



V« (J (.r./)f,C, 

,i,3 



-<c 



(2-13) 



since the expression in the middle of (|2.13|) is homogeneous of degree 
and the numerator and denominator are both bounded above and below 
on |£| = 1. From (|2.12|) we then have 



dxi 



dt 



< C 



which implies 



dx 



< Cv 7 ^. 



Integrating this inequality gives (1). 
2. Differentiating q with respect to x we have, 



J2( a tj)x(x,t)^j 

ij 



2 T,ai,j(x,t)&£j 



<C\£\ 



(2.14) 
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for some C independent of f , x, and t, since the expression in the middle 
of (|2 . 14[) is homogeneous of degree 1 and the numerator and denominator 
are both bounded above and below on |f | = 1 Using Gronwall's inequality 
gives 

— \— < |f (0)| exp(-Ci) < |f (i)| < |£(0)| cxp(Ci) < C(T, a) (2.15) 
6 (T, a) 

This results in the desired conclusion for hnitc T, that is for all t € [— T, T] 
if i < |f (0)| < a then there exists C(T, a) such that condition (|2.15l) holds. 

□ 

Recall that U(t,t') is the evolution operator associated to (|2.12l) and 

U(t,t )(x^y,f^y) (X^y(t, t , X^y, f^y), f^y(t, t ,X^y,fy)). 

By homogeneity, if c is a constant then the above equation scales as follows 

(Xry(t : t , Xj , <J%y) , C<J%y (t, £ , 3^ , f 7 ) ) — {x^yit^t , X 7 , Cf 7 ) , f 7 (t , t , 3^ , Cf 7 ) ) . (2.16) 

Since f 7 = 2 k uji^ with i < lu^l < 1 the relationship (|2.16p with c = 2~ k 
gives that the pair (x 7 ,w 7 ) lies in a compact subset of R" x (R"/{0}) whenever 
|x 7 | < R with R a constant independent of 7. We note that Lemma [4] then 
applies to (x 7 , w 7 ) so we that have a bound on the size of a; 7 (i) and f 7 = 2 k uji^{t) 
in terms of the initial data. 

Because our frame is similar to an almost orthogonal frame in type, it makes 
sense that the pairs of initial data which are close together in frequency con- 
tribute the most to the absolute value of the inner products in the sums in (|2.3[) 
and (|2.4|) . However, we have an extra variable a since we have a non-compactly 
supported set of frame functions. Therefore we will use the term 'close in fre- 
quency' to mean that the pairs (x 7 ,f 7 ) and (x 7 ',f 7 ') from Section 1 satisfy not 
only the condition |x 7 |, \x 1 / \ < R but also that \k — k'\ < fco, where fco is a finite 
constant independent of 7, 7'. We will show that close pairs of lattice variables 
have an extra property (Lemma [5j beyond that of Lemma [4] which makes it 
possible to compute the bounds on (|2.3[) and (|2.4|) . 

First we see by equation (|2.16[) . for all such close pairs with c = 2~ k 
(where here without loss of generality we have taken k' < k) the correspond- 
ing scaled pairs (x 7 , 2 k ~ k w 7 ) and (x 7 ',w 7 ') lie in the same compact subset 
[-R,R} n x [^2 k °] n of R" x (M"/{0». Thus we can conclude from Lemma H 
that the transformation U(t,0) is invertible and Lipschitz with uniform Lips- 
chitz constant when acting on (x 7 , 2 k ~ k w 7 ) and (a; 7 ',cjy). In other words, for 
all close pairs and for all t € [— T, T] with T fixed and finite there exists nonzero 
constants D\ and D2 independent of 7,7' with 

Z?id 2 ((x 7 (0,2' £ -' £ 'w 7 (t)); (ay(t),uy(t))) < d 2 ((x 7 , 2 fe " fe 'a; 7 ); 

(2.17) 

<D 2 d?((x y (t),2 k - k 'uj y (t)y,(x Y (t),u Y (t))) 
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where d denotes the usual Euclidean distance. For convenience, we will abbre- 
viate this type of equivalence relationship where the left hand side is bounded 
above and below by multiples of the right hand side by ~, so that inequality 
(|2.17|) can be rewritten as 

d 2 ((x 7 ,2 fc ~ fc 'u; 7 ); - d 2 ((x 7 {t), 2 k ~ k ' uo 1 {t)); (xy (t), ujy (t))). 

The inequality (|2.17|) allows us to obtain another, similar relationship which 
is the crucial relationship in the computations to obtain bounds on the action 
of the matrices and Bx(t). 

Lemma 5. For pairs (x 7 ,£ 7 ) and (a; 7 ',£ 7 ') such that |x 7 |,|x 7 '| < R where 
< R < oo and \k — k'\ < fco with R and fco independent of 7 and 7' the 
following holds 

d 2 (U(t,0)(x Y ,Lo y y i (x 7 ,2 fc - fc 'w 7 )) - d 2 (W(0,t)(x 7 ,2 fc - fc 'w 7 ); (x y ,lu y )) 
Proof. Since U(t, 0) o U(0, t) = I, the right hand side of the relationship, 
d 2 (W(0,t)(.T 7 ,2 fc - fc 'w 7 ); 

can be expressed as 

d 2 (U(0,t)(x 7 ,2 k - k 'uJ 7 );U(0,t)U(t,0)(xy,uJ Y )) 

and from estimate (|2.17p we obtain the desired conclusion. □ 

With these Lemmas, we can now calculate a bound on 

£M7,7',i)l (2.18) 
7' 

for fixed 7. We break this sum into three pieces: in region 1. 7' : k' < k — fco, 
in region 2, 7' : \k' — fc| < fco and in region 3, 7' : fc' > fc + fco where for all 
t G [-T,T],T < 00, fc = max{21og 2 C(T, a), 1}. For the rest of this argument 
let D > denote a constant which is independent of fc' and fc and is uniform 
for all t G [— T, T]. 

Wc will apply Lemma|4]in each region to subsets of the initial data (icy, wy) 
as outlined earlier. Here we must cutoff the a:y's so that |x 7 /| < i?, for some 
large positive R. This corresponds to having the initial data with support living 
in a ball of radius R. 

Again because of the similarity of the frame to an almost orthogonal frame, 
in regions 1 and 3 from Lemma 0] the exponential term from the bound on each 
inner product will dominate the sum, but in region 2 the argument is more 
subtle. In each case, formula (|2.16|) and Lemma |4] imply that: 

_ f \^ Y (t)\Ax 7 Ax Y y ( C?2S-< k 2%-' k '\^\\u> Y (t)\ \ f 
P7 ' 7 ' I Mi£yl + I£y'l) J \ 4^K| + 2*'K(0I) ) ^ > 



< 



C 2 27- £fe 2^- efc 'C(T,a) 
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The right hand side of (|2.1ip contains a product of two exponentials with argu- 
ments 



|2 fe w 7 - 2 fe 'w 7 /(i)| 2 

and 

l£yllM*)l ,„ - mi2 



(2.20) 



|£rl + l£r'(*)l 

In region 1, Lemma 0] implies a lower bound on (|2.19l) 



\x 7 - x Y (t)\ 2 (2.21) 



|2X-2 fc '^ft)| 2 (2 fc - 1 -2 fc 'C(T,q)) 



2 



4(2 fc |w T |+2 fc >y(t)|) 4(2 fc + 2 fc 'C(T,a)) 

2 fc - 5 (l-2 fc, - fc+1 C(T, a )) 2 2^ (l-c(fey) 2 fc 
> (l + 2«='-*=C(T, a )) > f 1 + _i\ ~ ~ 2 D - 

For (|2.2ip . we only know that 

\x n - Xy(t)\ 2 > 

which gives 

CXp (— \Xy — X 7 '(t)| 2 ) < 1. 

Since by assumption |a?y | < i? and jcy = Aa; 7 'a' = C e 2 _ T -efe <x' for fixed fc', by 

scaling we have \a'\ < RC e 2 ! ^' +ek . Bounding the number of points in 7L n in this 

ball by D(2 1 r +ek ') n , we obtain a bound on the number of xy for fixed (i, fc'). 
While the position of the x 7 < may change their total number does not change 

when they are propagated. From Lemma [T] there are 0{2~) vectors uii.k' in 
each annulus indexed by fc'. In region 1, from (|2.19|) we find since fc' < fc — fco 



C 2 2i -ek 2 ^--ek c{T,a) 



efea4- — efe' 



Combining estimates 



^ MT.V.t^O E 2-^-^2^+^ exp(-D2 



7':fc' <k~k Q I (z,fc') 

\fc <fe— feo 



(2.22) 

9 j ^ 2 nfc ' exp (-D2 fe ) J = O (2" fe exp (-D2 k )) 

\k'<k-ko / 
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But since 2 nk exp (— D2 k ) -> as k -> oo, the contribution from (|2.22|l is 
bounded independent of 7,7'. 

Similarly, in region 3 an application of Lemma 0] to the first part of the 
exponential contribution (|2.20p gives: 



|2 fc w 7 (0) - 2 k 'ujy(t)\ 2 (cp>)' 2fc ) 



2 



4(2 fc |w 7 (0)| + 2 k '\uj Y (t)\) A(2 k + 2 k 'C(T,a)) 

<-)k — 3 cyk— 2 

> > r>o k ' 

C(7» 3 C(T,a) 2 

Again, the same estimates as in region 1 for the number of the £y (t) and their 
exponential contribution hold, and the number of vectors co^k' in each annulus 

for fixed k' is still 0(2 s z~). For the size of /3 7 ,y from (|2.19p and the fact 
k' > k + ko, we have 



Thus 



/ 



£ |6 B ( 7 , 7 ',t)| = 

7':fc'>fe+fc 



£ 2T-^ 2 ^+^cxp(- J D2 A 



(2.23) 



= Ol J2 2 a ^+ i f+ 1 ^exp(-^2 fc ' N 

\fc'>fc + fe ' / 

By hypothesis, k' > k + ko, so the exponential term dominates the sum here as 
well so the contribution from (|2.23[) is bounded independent of 7,7'. 

If we try to simply apply Lemma 0] in region 2, as we did in regions 1 and 
3 we get a constant bound on the exponential contributions (|2.20[) and (|2.21[) 
which is not enough to dominate the contributions to the sum from the number 
of x 7 and £ 7 . Therefore the application of Lemma [5] to the exponential term 
is essential in region 2 since the treatment of the exponential contribution to 
the summation is more delicate there. The key is that the additional Lemma [5] 
allows us to sum over unpropagated variables which we have information about 
given from the construction that they are fixed in space. 

In region 2, by homogeneity and the fact |fc — k'\ < fco the entire exponential 
term can be re-written as follows 

M*)-er! 9 , r-^).^ (2 . 24) 



4(ierl + |£r'(*)l) |£rl + IM*)I 
-2 fc 'd 2 (W(t,0)(x 7 ,2 fc " fc 'a; 7 ); (x Y ,uJ Y )) 
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Applying Lemma (|5|) to (j2.24[) we obtain 

M*) "CI 2 \Cf{t)\\Cr\ , s r t) x |2 (2 25) 

4(ierl + l6y(*)l) + l6rl + IM*)l l ^ (J 7 ' { } 

where the constants in this equivalence relation may depend on ko as well but 
are uniform in T . 

Since xy = Ax 7 'a' we can factor out the scaling Axy = C e 2~~~ ek from 
part of the right hand side of (|2.25|) 

2 k ' \xy - x 7 (-0| 2 = C e 2 2" 2efc ' \a' - e | 2 

where we have set eo = (Axy')~ l x y (— t). Substituting A = C~ 2 2 2ek , an appli- 
cation of the integral estimates in the appendix gives 

]T cxp (-2 k ' \x Y - x 7 (-t)| 2 ) < CX P (-C 2 2~ 2ek ' \a' - e | 2 

a'6Z",|a'|<i? a'eZ™ 

(2.26) 

= o( 2 ek ' [ N 



But looking at inequality (|2.19[) the k! dependence in this last bound is exactly 
canceled by the size of /3 7 , 7 ' in region 2. Since \k — k'\ < kg, the other part 
of the exponential contribution we tackle with an argument similar to that of 
Lemma [1] applied to £ = 2 fe u; 7 (— t) which implies the sum 



^ exp I 2- k ' 2 k 'u y -2 k u y (-t) I (2.27) 



(i,k') 

\k'~k\<k 

is bounded independent of 7,7'. From these bounds and from the equivalence 
relation (|2.25p we can conclude 

M7,y,*)l = 0(l) (2.28) 

^y' :\k' — k\<ko 

and thus 

x;iM7,v,*)i=o(i) 
7' 

If we reverse the roles of 7 and 7' we can run a similar argument to the one 
above to bound 

£|M7,Y,«)I 
7 

The bounds in each of the regions \k — k'\ < ko, k > k' + ko and k < k' — ko will 
follow almost identically. The main difference in the argument will be the region 
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where \k — k'\ < ko we do not need to apply Lemma [5] since the 7 variables are 
not propagated. In this way we obtain the desired bound (|2.3|) . 
For the estimate (|2.4[) . we examine 



(2.29) 



The only difference between the bound on |&e(7, 7', t)\ and that of |6t(7, 7\t)\ 
is the factor of |£ n 
of Lemma 0] gives 



is the factor of |£ 7 — £ 7 '(£)| 2 + ICyCOII^yll 3 ^ — x ~t'{t)\ 2 - I R region 1, application 



l£r ~ £t Wl 2 + IM^II^IK - MO I 2 

< (lf 7 l + IM*)!)" + 1^7' WII^KK - 3* I + (*) - ^y I) 2 



< 



2 k + 2 k 'C(T, a) 



2 k 2 k' C(T,a)(R + T) 2 < D2 



2 k 



The rest of the estimates on /3 7j7 ' and the exponential contribution stay the 
same. Therefore by (|2.22j) 



E W7>y,i)i=o 

'j'-.k' <k — kg 



\ 



^2 2 2k 2 ~^-^ 2 ^+^ exp (_ D2 



■ (i,fe') 

\k'<k-k 



O I E 2 2fc 2 nfc ' exp (-L>2 fe ) j 

W<*s-&o / 



= O ( 2(™+ 2 ) fc exp (-D2 1 



and as k -> 00, 2 ( - n+2 ^ k cxp(-L>2 fc ) -> so the sum in question is also uniformly 
bounded independent of 7, 7' 

Similarly, in region 3, by Lemma 0] the extra factor is bounded by: 

|£y - £/ Wl 2 + M*) 1 1^7 I K " MOl' 
^1 + 1^7 



< (l£yl + l^7'(*)l) a + IMtJIfoKK - *7'l + IM*) ~ *y |) 2 



< 



2 fe + 2 k C ^ ^ + 2 k 2 k C( rp^ a ^ R + T ^2 < m 2k 

Again, the rest of the estimates stay the same, so that analogous to (|2.23j) 



E \b T (jn',t)\=0 

j':k'>k+k 



j2 2 2fe '2T-^2^+^exp(-^2 A 



/ (i.fc') 
■.k >fc — f-fco 



(2.30) 



0[ E 2 2fe W 

\fc'>fc + fc 



3nt' fcn 1 en(h'-h) 

• A • 9 



exp 



(-D2 fe ') 
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and as before this sum also converges independent of 7, 7' since k' > k + Icq. 

The only region where the extra factor in question makes a difference is in 
region 2. As in the treatment of the sum of |&e(7, 7')! over 7', Lemma [5] is 
again crucial. By homogeneity and Lemma [5] the extra factor in the bounds for 
|^r(7,7', t)\ can be rewritten as 



|e 7 -e7K0| 2 + I^WIIC7ll^-^Ki)| 2 ^2 2fe 'd 2 (W'(0,t)(.T 7 ,2^ fc ^ 7 );(xy, Wy )) 

and exponential factor in the bounds still follows the equivalence relation (j2.25[) . 
With these bounds in mind we split the sum 



\k—k' , 



E #y.7'(I£7 ~ £r'(*)| 2 + ier'(*)H6rll*7 - X T (*)| 2 ) 

7': I k—k' I <fco 



x cxp 



|2^ 7 -2*'av(t)| 2 |£ 7 ||M*)I 



ay(t)| 5 



4(2*|w 7 i + 2*'| Wy (t)i) ieyi + i^(*)r 

into two equivalent pieces. The sum (|2.31l) becomes 

>2fc' U . - t ,M 2 „™ (_ 2 fc ' \xy - X 7 {-t)\ 2 



E * 

7' :|fc— fc' I 



(— i)| exp 



x exp 2 



.-fc' 



2 fe 'u; 7 ' - 2 k LJ 1 (-t) 



and 



(2.31) 



(2.32) 



^2 Pin' 2 k 'uy - 2 fc w 7 (-i) 

7':|fe-fc'|<fc 

2 



exp ( —2 \x n 



x exp 2 



-A-' 



2 fe - 2 fe w 7 (-t) 



(2.33) 



To handle the sum (|2.32p , again since x 7 < = Ax 7 < a' we can factor out the scaling 



Ax Y = C e 2"- 



-ek' 



from part of the right hand side of (|2.25[) 



2 k ' \x. 



f (-t)\ 2 = C?2^ k ' \a' - e a \ 



and also from the new multiplying factor in (j2.32[) 

C 2 2 



2 2fe ' \x Y -x 7 {-t)\ 



2 _ r 2r ) ~2ek'+k' \ a , _ tQ f 



where we have set eo = (Axyi) 1 x 1 (—t). Substituing A = C e 2 2 2efc , an appli- 
cation of the second integral estimate in Appendix 2 gives 



( 22fe ' 1^7' - z 7 H)| 2 ) cxp (-2 k ' \xy - x^-t)\ 



a'& n ,\a'\<R 



(2.34) 



< 



E 



,-2ek'+k' il 



a - e \ exp -2 



-2ek' 1 / 



a - e 



^2*'(2 2efe ')t) 
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Using the previous estimates (|2.19p and (|2.27[) and the fact \k — k'\ < ko the 
sum (|2.32|) is 0{2 k ). For the second sum (|2.33|) , estimate (|2.26[) still applies for 
the sum over a' so we are reduced to examining 



|2 fc Wy -2XH) 



exp 



(i,k') 

| h — A" | ko 



2 k w 7 / - 2 fe w 7 (-i) 



(2.35) 



Now if we consider the same sets defined in Lemma Q] with £ = 2 k u!~ f (— t) for 
the first set A we get 



2 fc V 



2 k uj 1 {-t) 



< 2" 



implying from previous bounds on the number of cj 7 < in A 



^|2 fc 'w r -2 fc cj 7 (-<) 



exp 2 



2 fc w 7 , - 2 fc w 7 (-t) 



< 3 



n+loA; 



< 2ij + C(T,a)2 k < D2 k 



In each of the sets i3. 



2 k 'uj Y - 2 k u 1 {-t) 



< f2 k 



and similarly from an argument in Lemma 1 we can deduce 



2 fe w 7 (-t) 



exp 2 



2 k w y - 2 fe w 7 (-t) 



< ^2™2 fc (j + l)" +2 exp 



< D2* 



Now it is easy to see the contribution there is only a small (or 0) contribution 
coming from the sets C, T>, and £ (since \k — k'\ < ko), and this contribution is 
uniformly bounded independent of 7,7'. From here it follows that the second 
sum (|2.33|) is 0(2 k ). Combining the estimates on above gives 

E \b T (l^\t)\<D2 k . 

f':\k'-k\<k 

Since the contribution from regions 1 and 3 was uniformly bounded independent 
of 7, 7', we find 



By symmetry we can use similar estimates to obtain the second bound in 
2.41 Again, the main difference will be that there is no need to apply Lemma [5] 
in region 2. With these estimates we conclude the theorem. 
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3 Construction of the Parametrix 

With the frame of functions established, we turn our attention to constructing 
an appropriate parametrix for the Cauchy problem 

Tu(t,x) = (<9 t 2 - A(t,x,d x ))u{t,x) = 
u(t,x)\ t =o = f(x) 
d t u(t,x)\ t =o = h(x) 

where f(x) and h(x) are functions in L 2 (M. n ). We will construct operators 
C(t,t') and S(t,t') out of families of functions which are related to the frame 
functions. This section will follow the work of Hart Smith in [5] very closely. 

As earlier U(t,t') denotes the evolution operator associated to 'Hr = r — q. 
Additionally we denote the evolution operator associated to the Hamiltonian 
H+ = t + q as V(t, t'). Setting 

U(t, iO(x 7 (0), £,(())) = (x+(t,t'),£+(t,t')) 

and 

v(t,OMo),f 7 (o)) = (x-(t,t'),z-(t,t')) 

Then accordingly 
cf>±(t,t',x) 

= ( ) 2 exp (i$(t , t') ■ (x - x±(t, t')) - \$(t, t')\\x - s*(t, t')f) 

cf>±(t,t',x) 



and we let 

&(t,t',x) 



q(t' , Xj, £»y) 

From these definitions we construct the following operators C(t,t') 
Tl C(t,t')Tl f = P°B c (t,t')P?f 

= $>e(7,7',*M7')<M^) 

7,7' 

and S(t,t'): 

n s(t,t')u f = p°Bs(t,t')p°f 

= ^2bs(l 7 j',t)c(^')(l) 1 (x) 

7,7' 

where 



M7,7 , ,*)=o / ^(x)Ut(t,t',x)+<j)~,(t,t , ,x)) 



dx 
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and 

bs(7,l',t) = -J <j>^{x) (n+ (t, t', x) - n~, (t, t', x)) dx 

denote the entries of the matrices Bc(t,t') and Bs(t,t') respectively. 

Theorem 4. TC(t,t') and TS(t,t') are bounded operators of order one and 
zero respectively with operator norms which are uniformly bounded on intervals 
where t — t' is finite. 

Proof. By construction the matrix entries Bc{t, t') and B$(t, t 1 ) consist of linear 
combinations of inner products of propagated frame functions with the frame, 
so the proof is an immediate extension of Theorem [3] in section 2. □ 

Theorem 5. 

C(t',t')~I d t e(t',t') = o 

and 

S{t',t') = d t S{t',t')~I 

to leading order. 

Proof. The first statement follows directly from the definition of II and the 
integral calculations in section 2. For the second, we compute (where q = 
q(t',x l7 ^)) 

((d t n+(t',t',y)-d t n-(t',t',y))\ = ((l-^) (-1+*))^) 



q 2 ) 

as on null bicharacteristics r = ±q by homogeneity 

q 2 ^y 2 k ' 

and the result follows. From the proceeding arguments u(t, x) such that 

u(t, x) = S(t, t')h{x) + C{t, t')h{x) 
is the desired parametrix solution to the Cauchy problem. □ 

Theorem 6. If -1 < m < 2 and f e H m+1 (R n ),h e H m (W l ) and F £ 
L 1 ([-T,T],;# m (R 71 )) then there exists aGe Z^Q-T, T], ; H m (R n )) such that 

i 

u(t, x) = C(t, O)f(x) + S(t, O)h(x) + J (S(t, s)G(s, x) ds 

o 



([-T,T],;H m iTSt")) - C{T) ^||/|Ih™+i(R") + INIiJ™(R") + li 1 ([-T,T],;_H'"'(R»)) > 
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solves the Cauchy problem 

Tu(t, x) = (<9 t 2 - A(t, x, d x ))u(t, x) = F(t, x) 
u(t,x)\ t =o = f{x) 
d t u(t,x)\ t= o = h(x) 

in the weak sense. If f and h are both identically zero and F is also zero for 
yt G [— T, T] then G and u will vanish as well. 

Proof. As per Smith, we will show the existence of such a G using Volterra 
iteration. Assuming G e L l ([-T, T], ; H m (M. n )), we let 

v(t,x) = J S(t,s)G(s,x)ds 
o 

Because S(t, t 1 ) and dtS(t, t') arc both strongly continuous operators and S(t, t) = 
we have v{t,x) is in C([—T, T]; H m+1 (M. n )) n C l ([-T, T}; H m {W 1 )) and also 

d t v(t 7 x) = J d t S(t 7 s)G(s 7 x) ds 
o 

so it follows 

v(0,x)=0 d t v(t,x)\ t =o=0 
Furthermore differentiating in the sense of distributions we obtain 

d?v(t, x) = G{t 7 x) + J d?S{t, s)G{s, x) ds 
o 

We can conclude u{t, x) of the form in Theorem [6] is a weak solution to the 
Cauchy problem if the following Volterra equation 

G(t, x) + J TS(t, s)G(s, x) ds = F(t, x) - T (S(t, O)f(x) + S{t, O)h(x)) (3.1) 
o 

holds. Equation (|3 . 1 [) can be solved by iteration since the operator norm of 
S(t, s) is uniformly bounded on finite intervals of time by Theorem [3] Setting 

oo 

G(t,x) =F{t,x) + J2Gn(t,x) (3.2) 

n=l 

with 

t si s„_i 

G n (t,x) =[[... I S(t, Sl )S( Sl ,s 2 )... 







S(s„-i, s n )F(s„, x) ds n ... dsi 
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then G(t, x) is a solution to the equation 

G(t,x)+ ( S(t,s)G(s,x)ds = F(t,x) 



As the series in (|3.2[) converges in T, T], ; H m (W 1 )) with norm bounded 

by exp(TC(T) ||F||), this finishes the Theorem. □ 

Appendix 1 

It is well known that 

/e x p( % .„e x p(-^)*=© ? «p(-g) (3.3, 

R" 

We will use this fact to help us evaluate Integrals of the form 

J (y + b) exp (iy • 77) exp (-cy 2 ) dy (3.4) 

and 

J \y + b\ 2 exp (iy • 77) exp (-cy 2 ) dy (3.5) 

Recall that for c a constant, 77, y £ M" and a a multi-index with n components 
i\ a \ v a J exp(-cy 2 )exp(ir]-y) dy = (-1) |q| / <9^(exp (-cy 2 )) exp (iry • y) dy 

(3.6) 

and also 

exp (-cy 2 ) = -2cycxp (-cy 2 ) 

d 2 

exp (-cy 2 ) = (-2c + Ac 2 y 2 ) exp (-cy 2 ) 



With these equalities in mind, (|3.4|) is equal 
1 



<9 y (exp (-cy 2 )) exp (iry ■ y) dy + b / (exp (-cy 2 )) exp (777 • y) dy 
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We can also expand and re-write the integral in p.5[) so it is equal 
4c 2 / y - c 



^ dy(cxp (-cy 2 )) cxp (it) ■ y) dy+- I c\(cxp (-cy 2 )) exp (it) ■ y) dy 



+ (b 2 + J (exp (-cy 2 )) exp (ir] ■ y) dy 



Using the integration by parts formula 13.61 is just 

if \ ( v 2 ibn 

4c H 4c 2 c " 2c 



ex P "77 + — + b + 71 ( 3J ) 



4 Appendix 2 

Given an integer valued function h(a), h(a) may be estimated by the Euler 
summation formula 



r m p, 

£ h(a)= h{x)dx + Y,^V-\x)\tZ\ 

a<a<b i i=l 



i? m (4.1) 



where Bj is the j th Bernoulli number and hP(x) denotes the j th derivative of 
h(x). The remainder R m is defined as 

R m = (-l) m+1 [ ^Ml h m (x)dx 



The notation {x} denotes the fractional part of x, and B m ({x}) denotes the 
m th Bernoulli polynomial. Formula (|4.ip is derived in Concrete Mathematics, 
cf Ref[l]. 

Fix eo £ K n and A £ R, with A > 1. We wish to use the Euler summation 
formula to estimate the sums: 

E^UL^T) (4.2, 



and 

- P (-1^1!) (4.3) 

in terms of the parameter A. Since the variables 0:1,012, ...o n are indexed inde- 
pendent of each other, we may re-write (|4.2p as 



30 



We apply the Euler summation formula with to = 2 to the sum in parenthesis 
on the right hand side of (|4.4[) so that h(x) = cxp ^— — ^i! V Let g(x) = 
exp (—x 2 ) then via change of variables x = vA(u + e 0l ) 



R2 = lfJ B *({x})h"{x)dx= ^=jB 2 ({V\(u + e 0z )})g"(u)du 



By properties of the Bernoulli numbers (again, cf Ref[l]): 



\B 2 ({V\(u + e 0z )})\ <B 2 = l 



Integrating by parts gives 



J g"{u)du < J (4m 2 + 2)e- u2 du = 6^ 



\R 2 \ < 



The second term on the right hand side in the Euler summation formula van- 
ishes: 

m 



3- 



As a result 



E exp ( -l^L ]< (2^)1 



A 



The second sum (|4.3[) can be re-written as 

|"i ~ e 0i | 2 / |aj-e 0i | 2 
exp 1 



EE 

i=i \aiez 



A 



A 



E 



cxp 



,/ ,/ 12 



A 



Here a' = (a>i,a 2 , ...a t , ...a„) and £q = (e 0l , e 02 , ...e 0i , ■■•e „) Applying the Euler 
summation formula to 



\ a i- e 0i\ 2 ( |ai-£0iP 
r exp ' 



E 



gives (|4.5p is also 0{\T\). This follows since 



(4.5) 



x 2 fx 2 " 
— exp I — — ) dx = V 2ir\ 



(4.6) 



The details are left to the reader. Therefore (|4.3[) is C((A)S) as well. 
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